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ABSTRACT 

The Millennium Galaxy Catalogue (MGC) is a 37.5 deg 2 , medium-deep, _B-band imag- 
ing survey along the celestial equator, taken with the Wide Field Camera on the Isaac 
Newton Telescope. The survey region is contained within the regions of both the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey 
Early Data Release (SDSS-EDR). The survey has a uniform isophotal detection limit 
of 26 mag arcsec -2 and it provides a robust, well-defined catalogue of stars and galaxies 
in the range 16 < -Bmgc < 24 mag. 

Here we describe the survey strategy, the photometric and astrometric calibration, 
source detection and analysis, and present the galaxy number counts that connect the 
bright and faint galaxy populations within a single survey. We argue that these counts 
represent the state of the art and use them to constrain the normalizations (<t>*) of 
a number of recent estimates of the local galaxy luminosity function. We find that 
the 2dFGRS, SDSS Commissioning Data (CD), ESO Slice Project, Century Survey, 
Durham/UKST, Mt Stromlo/APM, SSRS2, and NOG luminosity functions require a 
revision of their published 0* values by factors of 1.05 ± 0.05, 0.76 ± 0.10, 1.02 ± 0.22, 
1.02 ± 0.16, 1.16 ± 0.28, 1.75 ± 0.37, 1.40 ± 0.26 and 1.01 ± 0.39, respectively. After 
renormalizing the galaxy luminosity functions we find a mean local bj luminosity 
density of = (1.986 ± 0.031) x 10 8 h L Q Mpc~ 3 .t 

Key words: catalogues - galaxies: general - galaxies: luminosity function, mass 
function - galaxies: statistics - cosmology: observations. 



1 INTRODUCTION 

Our understanding of the local universe and the local 
galaxy population originates primarily from the all- 
sky photographic Schmidt surveys and the established 
catalogues of bright galaxies derived from them, such 
as the Catalogue o f Galaxies and Clusters of Galaxies 
iZwickv et al. | Il96&f), the Morphological Cata lo gue o f 
Galaxies IVorontsov-VerYaminov fc Arkhipoval ta}74|h 
the Uppsala General Catalogue of Galaxies ([Nilson 

e , the ESO/Uppsala Catalogue jLaubertsI Il982t 
grt^_fe_Valerrtijn| LL9_89J), the Southern Galax y Cata- 
logue llCor^ri^leVaucouleurs fc de Vaucouleursl Il985t) . 



* E-mail: jol@roe.ac.uk 

f We use Hq = 100 h km s -1 Mpc - 1 throughout this paper. 



the Catalogue of Principal Galaxies iPaturel et al.lll989f) . 
the E dinburg h /Durham Southern Galax y Cat alogue 
jHevdon-Dumbleton. Collins fc MacGillivravl Il989l) . the 
APM catalogue jMaddox et al] ll990Jl. the Third Refer- 
ence Catalogue of Bright Galaxies ide Vaucouleurs et al, 
I l99lft and the SuperCOSMOS Sky Survey jHamblv et al. 
2001). While these catalogues have provided invaluable 
information and insight, uncertainty remains as to their 
completeness, parti cularly for low surface brightne s s and 
compact g alaxies ( Disney Il976t ISpravb errv et al] 119971 : 
llmpev fc BothuiJl997HDrinkwater et alJll999lh In addition 
there are concerns as to the photometric accuracy (e.g. 
lMetcalf e^Jbng^^_SJianl3LL995|), the susceptibility to scale 
errors llBertiii fc DennefeldliggTh . plate-to-plate variations 
llCross et al]l2003i) and dynamic range. 

These photographic-based catalogues have been the 
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starting point for numerous spectroscopic surveys aimed at 
measuring the local space density of galaxies (i.e. the local 
galaxy luminosity function). The space density of galaxies 
is our fundamental census of the local contents of space and 
there fore a crucial constraint for models of gal a xy formation 
(e.g. IWhite fc Frenkl Il99lt ICole et alJ l200d : iPearce et af] 
2001). If the imaging catalogues are in omission and/or pho- 
tometrically inaccurate then regardless of the completeness 
of the spectroscopic surveys our insight into the galaxy pop- 
ulation will be incomplete and most likely biased against 
specific galaxy types. 

Over the past two decades there have been nu- 
merous estimates of the local galaxy luminosity funo 
tion (e.g. EEP, lEfstathiou. Ellis fc Petersonl Il988l Mt 



Strom lo /APM iLovedav et all Il992t Autofi b. lEllis et al 
19961: ESP, IZucca et al] Il997t SSRS2. IMarzke et al 



199a Durham /UK ST. iRatcliffe et~ai] 1199 



, SDSS-CD, 

Blanton et al] 1200 it 2dFGRS. iNorberg et alJl2002ft and of 



_pa 

iSchechterl Il976h. . Typically the surveys agree broadly on 
the faint end slope (a, Aa « ±0.15) but show a marked 
variation in the characteristic luminosity (L* , AL* ~ 40 
per cent) and normalization ((/>*, Acj>* w 50 per cent). The 
uncertainties in the Schechter parameters result in an un- 
certainty of > 60 per cent in the local luminosity density, 

This uncertainty is usually expressed as the normaliza- 
tion problem which has been somewhat over shadowed by 
the more notoriou s faint blue galaxy problem llKoo fc Kronl 
Il992t lEllis! Il997ft . The latter describes the inability of 
basic galaxy number count models to predict the num- 
bers of galaxies seen at faint magnitudes (22 < B < 
28 mag) in the deep pencil beam CCD-b ased surveys 
(e.g. iTvsonl Il988t iMetcalfe et all Il995l 1200 lfl . The lesser 
known normalization problem describes the inability of 
number count models to explain the galaxy counts even 
at bright magnitudes (18 < B < 2 mag) by as much 
as a factor of 2 (see discus s ions i n IShanks et alJ ^84 
Driver. Windhorst fc Griffithsl Il995t IMarzke et aTT ll99S 



Cohen et aJl2003ft . In many ways the normalization problem 



is the more fundamental: while luminosity evolution, cosmol- 
ogy and/or dwarf galaxies can be, and have been, invoked 
in v arying mixtures to explain the faint blue galaxy problem 
(e.g. iBroadhurst. Ellis fc ShankJll988t iBabul fc Reedll992t 
iPhillipps fc DriveJll995t iFerguson fc Babullll99cJ) . none of 
these can be used to resolve the normalization problem. 

In the past the problem was typically circum- 
vented by renormalizing the number count mod els to 

the range 18 < B < 20 mag (e.g. |Djjver_et_al 

Metcalfe etafl Il995t [Driver. Windhorst fc Griffithsl 



Driver et al.ll99SUMarzke et alll998UMetcalfe et al 




The justification was that the bright galaxy catalogues, on 
which the luminosity function measurements are based, are 
shallow and therefore susceptible to local clustering. How- 
ever the crucial normalization range typically occurs at the 
faint limit of the photographic surveys (where the pho- 
tometry and completeness are more likely to be a prob- 
lem) and at the bright end of the pencil beam CCD sur- 
veys (where statistics are poor). While convenient, the clus- 
tering explanation overlooks two more worrisome possibili- 
ties: gross photometric errors and/or gross incompleteness 
in the local catalogues. If either of these two latter explana- 
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Figure 1. The magnitude ranges and survey areas spanned 
by some previous number count publications. Surveys based 
on photographic material are shown with a dashed line. 
The vertical lines show various transition regions where var- 
ious effects st a rt to dominate the galaxy c o unts. K ey: SSS 
iHamblv et alJ l200lft APM iMaddox etafl Il990fj|), SDSS- 
EDR iYasuda et all 1200 J). MAMA iBertin fc Dennefeldl 
Il997l) . EDSG C iHevdon-Dumbleton et all ll98Sft. MGC (thi s 



1991) 



work - ). G96 jGardn er et alJ Il996ft . J91 j Jones et alj 
EIS iPrandoni et alJ ll99Sft KW01 iKummel fc WagneJ 
186 Jlnfante. Pritchet fc Quintanal Il986t). K86 iKod .. 
A97 iArnouts et all Il997l)._ CADjS iHuang et al.l boOll). M91 



iMetcalfe et all 199 
l200ll) . M95 iMetcalfe ■ 1 



2001) 



1986) 



T88 jTvsonll988h , WHDF jMetcalfe et al 



HDF (Williams et alj|l99i 



al. 

96). 



1995ft . NTTDF iArnouts et*al 



199 



tions play a part this will have important consequences for 
the n ew-generation spect roscopic surveys , namely the 2dF - 
GRS dColless et al.ll200lTl and the SDSS llYork et alJl2000ft. . 
The input catalogue of the 2dFGRS is an extensively re- 
vised version of the photographic APM survey (which is 
known to show a pe culiar steepening in i ts galaxy counts 
at bright magnitudes. lMaddox et al.ll990bT) . with zero-point 
and scale-error corrections from a variety of sources includ- 
ing the 2MAS S Jf-band survey an d the data presented in 
this paper fsee lNorberg et al]l2002l for details). In the case 
of the SDSS - which leaps forward in terms of dynamic 
range, uniformity and wavelength coverage - the effective 
exposure time is relatively short (54 s) and the isophotal 
detection limit is comparable to that of the photographic 
surveys. Hence while issues of photometric accuracy should 
be resolved the question mark of completeness may remain. 

To address the above problems within a single, well- 
defined dataset we require a survey that is reasonably deep 
and yet has a large enough solid angle to provide accurate 
and statistically significant galaxy counts over the crucial 
normalization range. Furthermore, the survey's photometry 
must be accurate and its completeness high, i.e. it must 
probe to low surface brightnesses. Fig.^shows a number of 
imaging surveys in terms of their sky coverage and magni- 
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Figure 2. Outline of MGC fields 1, 2 and 3 and the arrangement 
of the CCDs. 



tude range. Dashed and solid lines indicate photographic 
and CCD-based surveys, respectively. Typically the faint 
surveys are CCD-based while the local surveys are photo- 
graph ic (with the notable an d recent exception of the SDSS- 
EDR. IStoughton et aljl2002ft . While the CCD surveys make 
significant improvements in surface brightness and magni- 
tude limits their sky coverage is small. It is only very re- 
cently that large CCD mosaics s uch as the Wide Field Cam- 
era (WFC.Ilrwin fe Lewisll200ll) and the SDSS instrument 
iGunn et~alT^)9gF) have been constructed that now allow a 
large area of sky to be surveyed within a realistic time frame. 

In this paper we present the Millennium Galaxy Cat- 
alogue (MGC, Sections |2H3- The MGC represents a new 
medium-deep, wide-angle galaxy resource, which firmly con- 
nects the local and distant universe within a single dataset 
(cf. Fig ID- in Section we produce the galaxy number 
counts spanning the range 16 < Bmgc < 24 mag. We then 
focus on the normalization problem by comparing our counts 
over the range 16 < -Bmgc < 20 mag to the predictions of 
a number of local luminosity function estimates in Section 
[3 Our counts provide stringent constraints on the normal- 
ization of the luminosity function and hence on the local 
luminosity density. Our conclusions are given in Section |H] 

The 2dFGRS and SDSS will essentially supersede all 
previous redshift surveys and therefore it is important to 
verify their photometric accuracy and completeness on as 
large a scale as possible. We will provide a detailed com- 
parison of the 2dFGRS and SDSS-ED R imaging catalog ues 
with the MGC in a companion paper IjCross et alJ l2003). 

A more long-term aim of the MGC project is to pro- 
vide structural information on the galaxy population around 
the crucial normalization point (16 < B < 20 mag). It is 
ironic that since the advent of the Hubble Space Telescope 
we have a greater understanding of the morphological mix 
of galaxies at faint magnit u des th an a t bright magn i tudes . 
For example, iDriver et all l|l99Sh and ICohen et all (|2003F l 
published morphological galaxy counts spanning the range 
21 < Bf45ow < 26 mag, yet no reliable morphological galaxy 
counts at brighter magnitudes exist. Consequently, no ac- 
curate local morphological luminosity functi ons exist (com- 
pare, for example, the confl icting results of lLovedav et all 
Il992l and lMarzke et aflll99Sl) and the evolution of the differ- 
ent morphological types cannot be accurately constrained. 
The MGC will enable us to remedy this situation as it allows 
morphological classification and the extraction of structural 
parameters to Bmgc = 20 mag. 
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Figure 3. Summary of the data quality across the MGC survey 
strip. The dots in the uppermost panel indicate the location of 
the photometric calibration fields. The sky and sky noise param- 
eters were calculated from the measured mode and rms of the 
background pixel value distribution. 



The data and catalogues presented in this paper are 
publically available at http://www.roe.ac.uk/~jol/mgc/ 



2 THE DATA 

2.1 The Wide Field Camera 

All data frames where taken using the Wide Field Cam- 
era (WFC). The WFC is mounted at prime focus on 
the 2.5 m Isaac Newton Telescope (INT) situated at La 
Palma. The WFC is a mosaic of four 4kx2k thinned 
EEV CCDs with a smaller 2kx2k Loral CCD which is 
used for auto-guiding. Each of the science CCDs mea- 
sures 2048 x 4100 pixels with a pixel scale of 0.333 arc- 
sec/pixel - this gives a total sky coverage of 0.287 deg 
per pointing. The four science chips are arranged as 
shown in Fig. |21 Full details of the WFC are provided 
at http://www.ast .cam.ac.uk/~wfcsur/technical.html (see 
also llrwin fe Lewisll200ll) . 

2.2 The observations 

The data constituting the MGC comprise 144 overlap- 
ping fields forming a 35 arcmin wide equatorial strip from 
9 h 58 m 28 s to 14 h 46 m 45 s (J2000). The observations were 
taken during 4 observing runs, 1999 March 15-16, 1999 April 
16-17, 1999 June 6-13 and 2000 March 26-April 4. Each 
field was observed for a single 750 s exposure through a Kitt 
Peak National Observatory B filter. Field 1 is centered on 
RA = 10 h 00 m 00 s , Dec = 00°00'00" (J2000) and field 144 
is centered on RA = 14 h 46 m 00 s , Dec = 00°00'00" (J2000). 
Hence each field is offset from the previous by 30 arcmin 
along the equatorial great circle. Fig. |5] shows the survey 
outline for the first three pointings. Note the substantial 
overlap between neighbouring fields. The survey region was 
chosen because it is contained within both the 2dFGRS 
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JColless et al.ll200ll) and the SDSS-EDR dStoughton et alJ 
2002) regions, thus providing redshifts and colours for the 
brighter galaxies and allowing a detailed check of the pho- 
tometry and completeness of these surveys. 

All observations were taken during dark or grey time 
through variable conditions. The seeing ranged from 0.9 to 
2.0 arcsec with the median seeing at 1.3 arcsec. The air 
masses rang ed from 1.141 to 1.672. Fig.Olshows a summary 
of the general observing conditions across the survey. 

As much of the data were collected during clear but 
non-photometric nights it was necessary to dedicate a sin- 
gle pristine photometric night (2000 March 30) to obtaining 
suitable calibration data at various stages along the survey 
strip. In total 20 MGC fields were observed during the pho- 
tometric night. Of these, six were only 100 s exposures as 



they had already been observed previously. These observa- 
tions were interspersed with 10 s observations of standard 
stars spanning a wide range in airma ss. The standard stars 
where taken from the lLandold (Il992l) standard areas SA98, 
SA101, SA104 and SA107. 

Some science frames were later found to be of too poor 
a quality to be useful and these were re-observed: the first 
eight fields were replaced with two 400 s exposures each and 
field 111 is a single 900 s exposure. 

2.3 Data reduction and astrometry 

All the preliminary data reduction - flat-fielding, bias cor- 
rection and astrometric calibration - was done by the Cam- 
bridge Astronomy Survey Unit (CASU) and full details of 
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Figure 5. (a) Median RA differences of doubly detected objects 
with 17 < B^GC — ^ mag in overlap regions. The error bars 
are 68 per cent ranges, (b) Frequency distribution of individual 
ARA values. This distribution has a mean and rms of —0.016 and 
0.083 arcsec, respectively, (c) Same as (a) for Dec differences, (d) 
Same as (b) for ADec. This distribution has a mean and rms of 
0.006 and 0.084 arcsec. 




1.2 1.3 


1.4 1.5 1.6 1.7 
sec Z 















0.2 0.4 0.6 0.8 1 1.2 1.4 

IB - V) 

Figure 6. Upper panel: colour-corrected zero-points versus air- 
mass. Lower panel: airmass-corrected zero-points versus colour. 
The solid lines show the best fit (cf . equation^ . The error bars in- 
clude contributions from photon counting, sky subtraction, read- 
out noise and the uncertainty in B. A further systematic error of 
0.01 mag (added in quadrature) was needed to achieve an accept- 
able fit. 



this process are provided bv llrwin fc Lewi's! (120011) . Briefly, 
a number of bias frames are collected each night and the 
median is subtracted from the data. All data (including 
flat-fields) are corrected for a known non-linearity. A twi- 
light fiat-field is taken during evening and morning twilight 
(when possible) and a median flat-field derived for that par- 
ticular run is divided into each data frame. After this process 
an initial astrometric calibration is made to the HST Guide 
Star Catalogue. Finally the frames are matched to the APM 
catalogue which itself is calibrated onto the Tycho-2 astro- 
metric system. Fig. 2] shows the final reduced image for one 
of our pointings, field 36. It illustrates problems with satel- 
lite trails, CCD defects, gaps between CCDs, etc. 

To assess the final astrometric accuracy we compared 
the positions of doubly detected objects in regions where 
neighbouring fields overlap (cf. Fig. The overlap regions 
are of size ~ 0.027 deg 2 and each contains ~ 60 objects in the 
range 17 < -B^ GC < 21 mag, 1 where detections can be easily 
and confidently matched. Fig.|S]shows the median positional 
differences for these objects for each overlap region as well as 
the overall ARA and ADec distributions. We find that both 
of these distributions have an rms of ±0.08 arcsec. Note, 
however, that they are slightly but significantly offset from 
zero which is most likely due to residual radial distortions. 



3 PHOTOMETRIC CALIBRATION 

As mentioned previously, four ILandohj lll992l) standard star 
fields were observed at a range of air masses throughout the 
course of the photometric night. For each observation of each 
standard star we computed a zero-point 



ZP S , 



B + 2.5 log/, 



(1) 



where B was taken from ILandohj <ll992l) and / is the flux 
of the star as measured from the data. We then fitted 
these zero-points with a double linear function in airmass 
(= secZ, where Z is the zenith distance), and colour, 

ZP std = a + a am sec Z + a B -v(B-V), (2) 

where (jB — V) is again taken from ILandohj £l992). In Fig.|S| 

we show the data and the fit. The residuals have an rms of 
0.02 mag and show no obvious trend with airmass, colour, 
B or time of observation. Note, however, that a systematic 
error of 0.01 mag was needed to achieve an acceptable fit. 
From the lower panel of Fig. HJ we can see that for several 
stars multiple observations of the same star give consistently 
high or low results. This may indicate that these stars are 
slightly variable or that the errors on the ILandohj <ll992l) 
photometry have been underestimated. 

For each MGC field we then computed a theoretical 
zero-point ZPth ~ a + a am sec Z, which is expected to be 
correct only for those fields observed during the photomet- 
ric night. We extracted objects as described in Section 2] 
and identified duplicate detections in regions where any two 
images overlapped. Using only objects with 17 < -Bmgc — 
21 mag and stellaricity > 0.5 (cf. Section l4.51 we computed, 
for each overlap region, the median of the magnitude differ- 
ences of the double detections, {A£?mgc)> as weu as an error 
on the median, ct(ab' mgc )- In Fig- a ) we show (AB^gc) 
for all the overlap regions using the initial, theoretical zero- 
points. 

A linear least-squares routine was then used to adjust 
the zero-points of the non-photometric fields in order to min- 
imize the quantity 



X 



E 



(ABmgc) 



(3) 



1 Here and in Section [3] BJ^GC re f ers to the KRON magnitude, 
uncorrected for Galactic extinction (cf. Section 14.31 . 



where the sum goes over all overlap regions. The process 
of object extraction, matching and zero-point adjustment 
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Figure 7. (a) Median magnitude differences of doubly detected 
objects with 17 < fi^GC — ^ mag and stellaricity > 0.5 in over- 
lap regions, using initial zero-points. Solid points indicate overlap 
regions involving at least one photometric observation. The error 
bars are 68 per cent ranges divided by the square root of the num- 
ber of objects (on average 35). (b) Frequency distribution of all 
of the individual AB' MGC values, (c) Same as (a) using the final 
zero-points. Note that the scale of the (ASj ll Q C )-axis is expanded 
by a factor of 10 compared with (a), (d) The solid line shows the 
frequency distribution of the individual AB' MGC values using the 
final zero-points. This distribution has a mean of —0.002 mag and 
an rms of 0.023 mag. For comparison, the shaded histogram is the 
same as that shown in (b). 



This is needed to reconcile the photometric fields 124 and 
139 with the photometric fields at field numbers < 74. 

Given the above calibration process the relationship be- 
tween an object's -Bmgc magnitude and its Landolt B mag- 
nitude is given by: 

B' MGC = B-0.U5(B -V), (5) 
where the error on the colour term is ±0.002. 



4 OBJECT EXTRACTION 

Object extraction was performed using Extractor, which 
is the S TARLINK adapted 2 versi on of S Extractor devel- 
oped bv lBertin fc Arnoutsl il996h . 



4.1 Background estimation 

S Extractor initially derives a background map by first 
defining a grid over the image and then passing a median 
filter (set to a size of 7 x 7 pixel) over each pixel within the 
grid cell. The local sky within the grid cell is then taken as 
the mode or cr-clipped mean of the pixel distribution within 
the cell. Finally, a background map is constructed via a bicu- 
bic spline interpolation over these points. We opted for the 
largest possible mesh size (256 x 256 pixel) to minimize the 
smoothing out of any extended low surface brightness fea- 
tures. 



was then repeated until a stable solution was reached (four 
iterations) . 

Since the zero-points of the photometric fields are held 
fixed at their theoretical values, the above procedure as- 
sumes that the observing conditions were perfectly stable 
throughout the photometric night. One can derive an esti- 
mate of the real-life error on the photometric zero-points, 
azp, by comparing the scatter of (AB MGC ) with cr( AS ; \ 
for those overlap regions that only involve photometric 
fields. We found azp = 0.005. Thus we modified the above 
calibration procedure by including the zero-points of the 
photometric fields in the parameters to be fitted and adding 
the additional constraints that they must lie 'close' to their 
theoretical values. In other words we now minimize the 
quantity 

all V ^ B MGO' / phot 

where the first sum again runs over all overlap regions and 
the second sum runs over all photometric fields. In Fig. fff c) 
we show (ABjj GC ) for all the overlap regions using the final 
zero-points. Fig. Eld) shows the histogram of the individual 
AB MGC values. The width of this distribution indicates an 
internal photometric accuracy of 0.023 mag for objects in 
the range 17 < Bmgc — 21 mag. Due to the paucity of 
photometric fields beyond MGC field 74 (cf. top panel of 
Fig. [3J the absolute calibration in the second half of the 
survey is less reliable than in the first: up to field 74 the 
median of the AB MGC distribution is 3 x 10~ 4 mag, which 
is consistent with zero, but beyond field 74 it is —0.004 mag. 



4.2 Object detection and deblending 

After convolving the image with a filter, SExtractor de- 
tects objects as groups of connected pixels above the uni- 
form surface brightness threshold of jUu m = 26 mag arcsec -2 . 
Once an object has been detected SExtractor redetects 
the object using 30 different detection thresholds, which are 
exponentially spaced between the peak flux value of the ob- 
ject and the detection threshold. If at any level the object 
breaks up into two or more disconnected subcomponents, 
each containing at least 10 per cent of the total flux of the 
object, then the object is deblended. This multithresholding 
is ideally suited for galaxy extraction because no assump- 
tions concerning the sh ape of the object are being made 
jBertin fc Arnoutslll996T) . 

4.3 Photometry 

The photometry was performed using SExtractor with a 
constant analysis isophote of 26 mag arcsec -2 to provide a 
uniformly processed catalogue. Four of the various types of 
magnitudes provided by SExtractor are included in the 
MGC: an isophotal magnitude (ISO), a corrected isopho- 
tal magnitude (ISOCOR), an adaptive aperture magnitude 
(KRON) and a best magnitude (BEST). The ISOCOR mag- 
nitude is calculated by correcting the ISO magnitude for the 

2 Note that the STARLINK version has additional data handling 
routines, the STARLINK World Coordinate System software and 
a graphical interface. In all other aspects the code is ident ical to 
that developed and described by iBertin fc Arnoutsl Jl996l) , thus 
henceforth we refer to SEXTRACTOR in the text. 
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Figure 8. Random examples of the field galaxy population at -Bmgc = 16.25 mag (left) in steps of 0.5 mag to Bmgc 
(right) displayed from 22 mag arcsec -2 (black) to 3<r below the sky (white). The image sizes are 33 X 33 arcsec 2 . 



19.75 mag 



fraction of the flux outs ide the limiting isophot e assuming 
a Gaussian profile (see iBertin fc Arnoutslll996l . and refer- 
ences therein for details). KRO N magnitud es (i.e. ellipti- 
cal apertures of 2.5 Kron radii, lKrorJll980t) are known to 
underestimate the fluxes in perfect exponential profiles by 
0.04 mag and in de Vaucouleur profiles by 0.1 mag. Never- 
theless, they have been shown to be the most robust to vari- 
ations in redshift, bulge-to-disc ratio, isophotal limit and 
seeing llCrossl 12002*) . The BEST magnitude is taken to be 
the KRON magnitude except in crowded regions where the 
ISOCOR magnitude is used instead. All objects are individ- 
ually corrected for Gal actic extinction using the dust e xtinc - 
tion maps provided bv lSchlegel. Finkbeiner fc Davisl dl99ftl) 
and adopting Ab kpno = 4.23. From now on we refer to the 
BEST magnitudes before and after extinction correction as 
Bmgc an d Bmgc, respectively. Fig. |^1 shows a random se- 
lection of galaxies in the range Bmgc = 16.25 to 19.75 mag. 



4.4 Overlap regions 

As a result of the substantial overlap regions (cf. Fig. the 
catalogue contains many duplicate objects. These were used 
in previous sections to verify the astrometry and to cali- 
brate the photometry. We now remove duplicate detections 
by imposing RA limits for each field, effectively splitting 
each overlap region in half. (Note that the RA limits do not 
apply to objects detected on CCD 3.) For a small number of 



objects, all lying very close to an RA limit, this procedure 
did not remove one of the duplicate detections, because the 
two detections happened to lie on either side of the limit. 
These cases were fixed by hand. 



4.5 Classification and cleaning 

At this stage the catalogue contains a total of 1 070 374 ob- 
jects to Bmgc = 24 mag. The catalogue comprises galaxies, 
stars and various unwanted objects and artefacts such as 
satellite trails, CCD defects, cosmic rays, diffraction spikes, 
asteroids and spurious noise detections. As a starting point 
for classification we used the stellaricity parameter pro- 
vided by SExtractor, which is produced for each object 
by an artificial neural network (ANN) that has been ex- 
tensively trained to diffe rentiate between stars and galaxies 
jBertin fc Arnoutslll996l) . The input of the ANN consists of 
nine object parameters (eight isophotal areas and the peak 
intensity) and the seeing. The output consists of a single 
number, called stellaricity, which takes a value of 1 for stars, 
for galaxies and intermediate values for more dubious ob- 
jects. Fig. |U] shows the number of objects as a function of 
stellaricity and B' MGC . At B^gc ^5 20 mag the stellaricity 
distribution is clearly bimodal with almost all values at the 
extremes and so star-galaxy separation is trivial. At fainter 
magnitudes the star-galaxy separation requires more effort. 
Hence we now define two catalogues: MGC-BRIGHT, which 
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Table 1. Breakdown of MGC-BRIGHT (B M GC < 20 mag). 




0.5 1 1.5 2 2.5 5 3.5 4 

Figure 9. Number of all detected objects as a function of stel- 
laricity and B^ GC . Dots indicate the positions of objects from 
the field with the very worst seeing. The vertical dashed line 
marks the approximate division between MGC-BRIGHT and 
MGC-FAINT (which was defined in B MGC , not B^gc)- In MGC- 
BRIGHT, the solid line delineates the division between objects 
assumed to be stars and those classified by eye. In MGC-FAINT 
the three solid lines show the star-galaxy separation adopted for 
the fields with the best, median and worst seeing (see the text for 
details), where the uppermost line corresponds to the best seeing. 



contains all objects with Bmgc < 20 mag, and MGC-FAINT 
which contains the rest. 



4.5.1 MGC-BRIGHT: B MG c < 20 

For MGC-BRIGHT we adopted the following classification 
strategy. Brighter than B[j GC = 15 mag we classified all 
objects with stellaricity > 0.65 as stars. This low stellar- 
icity cut was used because the objects in the upper left- 
hand corner of Fig. |5] are, in fact, flooded stars. For the 
rest of MGC-BRIGHT we classified all objects with stel- 
laricity > 0.98 as stars. This is a 'natural' value to adopt 
because the stellaricity distribution rises sharply from 0.97 
to 0.98. All objects so far classified as non-stellar were then 
inspected visually and classified into one of the following 
categories: galaxy, star, asteroid, satellite trail, cosmic ray, 
CCD defect, diffraction spike or spurious noise detection. 
Incorrectly deblended galaxies were then repaired by hand, 
but their original catalogue entries were retained and clas- 
sified as obsolete. Asteroids were verifi ed using images from 
the SuperCOSMOS Sky Survey (SSS. lHamblv et aljEoOll) . 
We also inspected all objects with FWHM, semimajor or 
semiminor axis less than the seeing (these turned out to be 
primarily cosmic rays and CCD defects). During this pro- 
cess galaxies were also assigned one of three quality classes, 
Q, depending on the level of: (i) contamination by CCD de- 
fects, satellite trails, cosmic rays and diffraction spikes; (ii) 
blending with a similarly bright object; (iii) missing light 
due to a CCD edge and (iv) failed background estimation 
due to nearby bright objects. The breakdown of the final 
MGC-BRIGHT catalogue is shown in Table Q and is a good 
indication of the level of contamination in purely automated 
galaxy catalogues. 



Description Number After cleaning 



Galaxies 11866 9913 

(Q a = l,2,3 11266,449,151 9775,138,0) 

Stars 51284 42365 

Asteroids 145 125 

Satellite trails 162 

Cosmic rays 113 62 

CCD defects 3027 

Diffraction spikes 263 

Noise detections 2023 13 

Obsolete 140 116 

Total 69023 52594 



"Quality class, where 1 denotes highest quality. 
4-5.2 Exclusion regions 

Having classified MGC-BRIGHT we now have reasonably 
good indicators for the positions of CCD defects, satellite 
trails, very bright objects and diffraction spikes. All of these 
adversely affect the measurement of parameters of nearby 
objects. In addition, very bright objects (Bmgc < 12.5 mag) 
cause a halo of spurious faint detections. Before we continue 
with the star-galaxy separation in MGC-FAINT it is impor- 
tant to remove as many of these adversely affected and spuri- 
ous detections as possible. Hence we now define exclusion re- 
gions: any objects within these regions will be removed from 
both MGC-BRIGHT and MGC-FAINT to produce 'cleaned' 
versions of these catalogues. 

In particular, we define rectangular or elliptical exclu- 
sion regions around CCD edges, the vignetted corner of CCD 
3, CCD defects, CCDs 3 and 4 of field 79 (which failed to 
read out properly) , small, unwanted overlaps between CCDs 
2 and 3 of a few neighbouring fields, satellite trails, diffrac- 
tion spikes and objects with £?mgc < 12.5 mag. For most 
'classes' of exclusion regions we use some simple algorithm to 
define a first set of exclusion regions from the correspond- 
ing class of detections in MGC-BRIGHT. This first set is 
then improved upon and augmented by hand where neces- 
sary. Since the parameters of very bright objects in MGC- 
BRIGHT are unreliable due to saturation, we have primarily 
used the SSS to define the last class of exclusion regions. 

Fig. |l]shows the exclusion regions for field 36. The ex- 
clusion regions reduce the total area of the survey from 
37.50 deg 2 to 30.84 deg 2 . 

Fields 14, 15 and 65 are of substandard quality be- 
cause their surface brightness detection limit is considerably 
brighter than 26 mag arcsec _2 (cf. bottom panel of Fig. |3J. 
This results in a large number of spurious detections at mag- 
nitudes > 23 mag in these fields. As part of the cleaning pro- 
cess we have therefore removed all objects from these fields 
from MGC-FAINT (but not from MGC-BRIGHT). 

4.5.3 MGC-FAINT: Bmgc > 20 

Although MGC-FAINT has been cleaned as described 
above, it still contains large numbers of cosmic rays, which 
we now attempt to identify. Cosmic rays are expected to be 
very small along at least one axis. SExtractor provides 
object shape parameters in the form of an 'rms ellipse'. In 
Fig. llOl we plot all objects from the cleaned versions of MGC- 
BRIGHT and MGC-FAINT in terms of their flux rms along 
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Figure 10. Flux rms along the minor axis versus B' MGG for all 
objects from the cleaned versions of MGC-BRIGHT and MGC- 
FAINT. Almost all cosmic rays lie in the band delineated by the 
two solid lines. For B' MGG < 22.7 mag the band is well separated 
from the general population and hence reliable and complete cos- 
mic ray identification is possible using a minor axis cut alone. 
At B' MQC > 22.7 mag we apply the additional cut of FWHM 
< 0.9 X seeing to identify cosmic rays but the selection is now 
incomplete. 



the minor axis and B MGC . We can clearly identify a 'band' 
of objects with 0.7 < minor axis < 1.0 pixel (solid lines). 
From our classification of MGC-BRIGHT and inspection of 
random samples at fainter magnitudes we found that all ob- 
jects within this band and B' MGC < 22.7 mag are indeed 
cosmic rays. Conversely, clearly almost all cosmic rays lie 
in this band. However, at B M gc ^ 22.7 mag the cosmic ray 
band merges with the general population of objects and a 
minor axis cut alone is insufficient to select cosmic rays. 

S Extractor provides an estimate of the FWHM of 
the flux profile, assuming a circular Gaussian core. The ad- 
ditional cut of FWHM < 0.9 x seeing is successful in identi- 
fying large numbers of cosmic rays at Bmgc > 22.7 mag. 
The remaining unidentified cosmic rays are those with a 
large angle of incidence which leave a faint 'trail' on the 
images. Unfortunately, given the parameters available from 
SExtractor this class of cosmic rays is genuinely indistin- 
guishable from real objects, and hence the identification of 
cosmic rays at B' MGC > 22.7 mag must remain incomplete. 
In Section 16.41 we estimate this incompleteness to be < 20 
per cent and correct the galaxy number counts accordingly. 

All objects in the cleaned version of MGC-FAINT not 
identified as cosmic rays by the above procedure are as- 
sumed to be either stars or galaxies. As we have noted 
above, for B M gc ^ 20 mag star -galaxy separation b e comes 
increasingly difficult. However, iKiimmel fc Wagner! j200ll) 
showed that the slope of the star counts remains constant 
to B — 22.75 mag, in agreement with Galax y model predic- 
tions feahcall fc Soneiral 19801 iBahcall 198r3l . Hence one can 



produce statistical galaxy number counts without explicitly 
classifying individual objects by simply subtracting the ex- 
pected stellar counts (determined by extrapolation from the 
bright end) from the measured total counts. 

Nevertheless, a variety of applications do require classi- 
fications for individual objects and so we adopt the following 
scheme. The idea is to find a stellaricity value, s c , such that 
the objects with stellaricity > s c reproduce the numbers of 
stars expected by extrapolating the stellar counts measured 
in the range 17 < Bmgc < 20 mag. However, the abil- 
ity of SExtractor to distinguish between stars and galax- 
ies depends both on magnitude and the seeing. Moving to 
fainter magnitudes or worse seeing both have the effect of 
redistributing stars from the high end of the stellaricity dis- 
tribution towards intermediate values. Therefore, we must 
determine s c as a function of Bmgc f° r each MGC field in- 
dividually. 

Since for a given field s c (B M gc) should be monotonic 
and since the exact placement of s c is most crucial at the 
faint end, we have adopted the following strategy: for each 
field we first adjust s c in the field's faintest three bins (de- 
termined by the field's completeness limit, see Section [4.61 
to give star count values closest to the extrapolated bright 
star count fit (in a x 2 sense), while requiring s c (B M gc) to 
be monotonic. In all remaining bins with B' MGC > 20 mag 
we set s c to the highest s c value so far obtained. The bright 
star count fit is derived locally for each field from the stel- 
lar counts in the range 17 < B' MGC < 20 mag as measured 
from the five fields centered on the field under considera- 
tion (to improve the reliability of the fit). In Fig.|U]we show 
the resulting s c curves for the fields with the best, median 
and worst seeing. As expected we found that the value of s c 
correlates strongly with seeing. 

The error on the galaxy number counts introduced by 
the uncertainty in the exact placement of s c is discussed in 
Section 16.11 



4.6 Faint completeness limits 

In order to estimate the faint completeness limit for each 
field we used the artdata package of IRAF to add 100 
equally bright point sources at random positions (but avoid- 
ing exclusion regions) to each CCD of each field. We then 
re-extracted object lists to determine whether the simulated 
objects were recovered. Detections that were judged to be 
due to a blend of a simulated object and a real object of 
similar or greater brightness were not counted as recovered. 
This process was repeated for a total of 14 input magni- 
tudes in the range 22 < B M gc < 26 mag. For each field we 
then fitted a three-parameter completeness function to the 
fraction of recovered obje cts as a function of magnitude (cf. 
IKiimmel fc WagnerlEoOlh : 



c (Bmgc) = 



exp 



Bmgc ~ B'so 



1 



(6) 



where Cb is the completeness at bright magnitudes and B' 50 is 
the magnitude at which the completeness reaches c(B5 ) = 
Cb(cb + ~ 0.5. In general ^ / 1 because some small 
fraction of objects will always be covered up by brighter 
objects and thus go undetected. 
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Figure 11. The distribution of galaxies and stars in MGC- 
B RIGHT (top panel) and MGC-FAINT (bottom panel) as a 
function of MGC field number or RA (indicated along top). The 
dashed vertical lines indicate the positions of Z < 1 galaxy clus- 
ters found in the literature. The galaxies show a uniform distri- 
bution whereas the stars show a clear rise towards the Galactic 
Centre. The solid horizontal lines indicate the mean number of 
bright and faint galaxies per field. The curved solid lines repre- 
sent the star counts pred icted by the standa rd Galaxy model of 
iBahcall fc Soneiral ll98Cft and iBahcall <1986t) . All numbers have 
been scaled to a constant area of 0.215 deg 2 per field. The dot- 
ted lines and right axes indicate how the extinction and galactic 
latitude vary across the survey. 



We find that Cb is weakly correlated with seeing, varying 
from 0.995 to 0.965 with a mean value of Cb = 0.983. B' 50 is 
quite strongly correlated with seeing while f3 correlates with 
the sky noise (f) = 0.22). 

For each field we now define a magnitude limit by 
c(B[ im ) = /cb, where we take / = 0.97 so that the mean 
completeness at the magnitude limit is c(B[ irn ) — fcb > 
0.95. Finally, for a given field we define a dust-corrected 
magnitude limit, -Bii m , as the brightest Bmgc of all objects 
with B' MGC > B' liul in that field. In the following we will only 
use objects with Bmgc < -Biim and we will not attempt to 
use fainter objects in combination with incompleteness cor- 
rections. Hence, when constructing the number counts we 
will only use a given field for a given magnitude bin if the 
field covers the entire bin (of size ABmgc = 0.5 mag). Thus 
we are effectively rounding the magnitude limit of the field 
downwards to the nearest multiple of 0.5. This leaves 124 
fields in the Bmgo = 23.25 mag bin and 14 fields in the 
faintest bin at Bmgc = 23.75 mag. 



4.7 Object distribution in RA 

Fig. [HI shows how the MGC-BRIGHT (upper panel) and 
MGC-FAINT (lower panel) galaxies are distributed along 
the MGC survey strip. The vertical dashed lines indicate the 
positions of known z < 1 galaxy clusters. It is encouraging 
to note that on the whole they coincide with peaks in the 
galaxy numbers. 

We also show how the numbers of stars vary across 
the survey strip and how they compare to predictions of 




Figure 12. The apparent bivariate brightness distribution of 
MGC-BRIGHT galaxies. The shaded regions show the selection 
boundaries determined by the median seeing (upper diagonal 
line), the magnitude limit (vertical line), the surface brightness 
detection limit (horizontal line) and the mesh size employed dur- 
ing the background estimation (lower diagonal line, see the text 
for details). The two dashed lines are the selection limits for the 
fields with the worst (lower) and the best (upper) seeing. 



the standard G alaxy model of lBahcall fc Soneiral lll980l) and 
IBahcall l)l986l) . The basic shape of the counts and the model 
agree in that there is a clear rise towards the galactic bulge. 
However, there appear to be significant differences. The 
model systematically under-predicts the faint counts and 
there are also discrepancies with the bright counts at low 
galactic lati tudes. Similar trends are apparent in the stellar 
r* counts of lYasuda et alJ l|200ll ) . We will address th is issue 
in more detail in a future paper (iLemon et al]l2003l) . 



5 SELECTION LIMITS 

One of the primary aims of the MGC was to define a sam- 
ple of galaxies with well-defined selection criteria, hence the 
use of a constant detection isophote. Fig. 1121 shows the dis- 
tribution of galaxies in the apparent magnitude - apparent 
effective surface brightness plane for MGC-BRIGHT. The 
apparent effective surface brightness was simply calculated 



/Ueff = Bmgc + 2.51og(27rr 1/ - 2 ) 



(7) 



where rx/2 is the semimajor axis of the ellipse containing 
half the total flux (in arcsec) which was measured directly 
from the data. 

The solid lines delineating the shaded regions in Fig. 1121 
show the selection boundaries. The three principal se- 
lection limits are the median seeing limit of 1.3 arcsec, 
the central surface brightness detection limit of /nu m = 
26 mag arcsec -2 (equivalent to an effective surface bright- 
ness limit of 27.12 mag arcsec -2 , assuming an exponen- 
tial profile) and the imposed magnitude limit of Bmgc = 
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20 mag. Note that this cut was set by the magnitude limit 
at which star-galaxy separation is reliable. 

A fourth selection limit is implied by the method of 
background estimation described in Section r4.ll Any objects 
covering a substantial fraction of a cell of the grid used to 
estimate the background will produce an erroneously high 
background measurement and hence may be missed. We es- 
timate that objects with an isophotal area of ~ 25 per cent 
of a grid cell will be affected. We have used a mesh size of 
256 x 256 pixel (85.2 x 85.2 arcsec). Assuming n so w 2ri/2 
this represents an upper size limit in terms of half-light ra- 
dius of 64 pixel (21.3 arcsec). This line is shown as the lower 
diagonal selection limit in Fig. 1121 

Outside the selection boundaries galaxies cannot, theo- 
retically, be detected although noise may scatter a few ob- 
jects across these boundaries. Galaxies with parameters in- 
side the boundaries should be detectable. The galaxy pop- 
ulation follows a well-defined distribution, which does not 
reach to the high and low surface brightness selection bound- 
aries, demonstrating the robustness of MGC-BRIGHT with 
respect to surface brightness selection effects. The distri- 
bution in surface brightness of the observed population is 
far too narrow to b e explained by visibility theory (see 
ICross fc Driverll2002T) and one must conclude that luminous 
low surface bright n ess galaxies a r e inde ed rar e as suggested 
bv iDriveJ ll99Sh . ICross et alJ (1200 J) and iBlanton et ail 
fcOOll) . 
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Figure 13. Various estimates of the error on the MGC galaxy 
number counts, ctjv, scaled by the number counts, -/VmgCi as i n ~ 
dicated and described in the text. The last two sets refer to the 
MGC-2dFGRS counts discussed in Section The crosses and 
'+' symbols show the final errors listed in Table l2j 



6 NUMBER COUNTS 
6.1 Errors 

Before we present the final galaxy counts we will estimate re- 
alistic errors for the counts. In addition to Poisson noise they 
should also contain a contribution from large-scale structure 
(LSS), which will also induce correlations between differ- 
ent bins. If galaxies have an angular correlation function 
w(0) — A W J 1 then the error on the counts in a given bin, 
crjv, is given by dPeeblesl lT980) 

a N (B M ac) = A^y/N + Af 2 A w C, (8) 

where J\f is the number of galaxies in the bin, A is the sur- 
vey area and C is the integral constraint given by a double 
integral over the survey area: 

C = A' 2 J J ^-Mwidwa, (9) 

where we will use the SDSS-EDR results on w(ff) of 
IConnoliv' et alJ J2002I) . who found 1 - 7 » -0.7 and mea- 
sured log A w as a function of r* . We translate these mea- 
sureme nts to -Bmgc (using equation the colour equa- 
tions of |Fukugit a et aljfl996l and the mean galaxy colours of 
lYasuda et alJl200ll) and extrapolate where necessary. The 
resulting error estimates are shown as cross es in Fi g. 1131 
In fac t, we follow the detailed description of lYasuda et alJ 
( 2001) to calculate the full covariance matrix for the counts, 
which will be needed in Section 17.31 

We have also e stimated the covariance matrix using the 
jackknife method (Efron 1982]). This was implemented by 
excluding each field from the dataset in turn and measuring 
the counts from the remaining data. The covariance of these 



Nf measurements is then multiplied by (JVf— 1) , where 
Nf is the number of fields used for a given pair of bins. The 
result agrees well with the integral constraint estimate above 
and we show the resulting errors as open squares in Fig. 1131 

Finally, we have estimated the errors using the 2dFGRS 
imaging catalogue (2dFGRS-IC). We have divided the North 
and South Galactic Pole regions (NGP and SGP) into 21 and 
32 subregions, respectively, each of similar size and shape as 
the MGC survey region. We converted the 2dFGRS-IC mag- 
nitudes to the Bmgc filter system using -Bmgc = bj + 0.074 
(which includes a known offset of 0.056 mag; cf. eauation ll3l . 
We then measured the counts in each independent subre- 
gion and plot the rms of these sets separately for the NGP 
and SGP as triangles in Fig. 1131 This is only possible at 
Bm gc < 19 niag due to the magnitude limit of the 2dFGRS- 
IC JColless et alj|200ll) . 

We note that all of the above error estimates are in 
reasonably good agreement and lie well above the Poisson 
errors (solid dots). Henceforth we will adopt the integral 
constraint estimates. 

In Fig. 1131 we also plot as open circles an estimate of 
the uncertainty introduced by the star-galaxy separation in 
MGC-FAINT. We have varied the method of determining 
s c (cf. Section [4.5.31 in numerous ways but we only found 
a significant change in the number counts by systematically 
increasing (or decreasing) s c in all fields. The error shown in 
Fig. 1131 is the difference in the counts arising from a global 
change of s c by 0.01. We note that such a change produces 
a severe mismatch between the measured faint star counts 
and those expected from extrapolating the bright counts. 
Hence the derived error is probably an overestimate. In any 
case, it is well below the jackknife errors in all but one bin 
and we therefore do not consider it any further. 
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Figure 14. The top panel shows various galaxy number count 
estimates derived from the 2dFGRS-IC as indicated: full sample, 
NGP only (740.3 deg 2 ), SGP only (1094.7 deg 2 ), in the MGC 
survey region, in 21 independent NGP subregions, each of a size 
and shape similar to the MGC region. We also plot the counts 
predicted by the 2dFGRS luminosity function (cf. Section \7. II . 
The error bars on the NGP and SGP counts were calculated using 
the integral constraint method. The middle panel shows the LSS 
correction factors using the full, NGP only and predicted counts 
(same symbols as in top panel). The bottom panel shows the 
uncorrected MGC counts (crosses) as well as the MGC-2dFGRS 
counts using the correction factors of the middle panel (same 
symbols). Several sets of counts in the top and bottom panels 
have been offset by ±0.1 mag for clarity. 



6.2 Correcting for large-scale structure 

Since the MGC survey region is contained within the 
2dFGRS-IC NGP region the question arises as to whether 
the photometric accuracy and high completeness of the for- 
mer can be combined with the large area of the latter to 
derive number counts less affected by LSS and hence of 
higher accuracy. The idea is to view the ratio of iVic-MGC> 
the 2dFGRS-IC counts within the MGC region, to Nic, the 
counts from a much larger 2dFGRS-IC sample, as an LSS 
correction factor, which can be applied to the MGC counts 
to give LSS-corrected counts: 



«mgc — 



IC-MGC 



NmGC- 



(10) 



Equivalently, one can view Nmgc/Nic-mgc as a factor be- 
ing applied to Nic, which corrects for incompleteness, stel- 
lar contamination, bad deblending, etc. in the 2dFGRS-IC. 
Hence we will refer to A^mgc as the MGC-2dFGRS counts. 

Fig. 1141 illustrates the LSS correction procedure. In the 
top panel we plot JVic.mgc as crosses along with four sets 
of large-area counts for possible use as Nic (open symbols). 
The middle panel shows the LSS correction factors using 
the full, NGP and predicted counts (same symbols as in top 
panel), and the bottom panel shows both the uncorrected 



MGC counts (crosses) and the MGC-2dFGRS counts using 
the above correction factors. 

For this correction procedure to be valid we require only 
that the properties of the 2dFGRS-IC (photometry, incom- 
pleteness, etc.) are homogeneous over the region used to 
derive Nic and the MGC region. However, in the top panel 
of Fig. I14l the SGP appears very significantly underdense in 
comparison to the NGP. Indeed. fNorberg et al] l)2002l) found 
a 7 per cent difference in the numbers of NGP and SGP 
galaxies to bj — 19.2 mag but concluded that this was 'rea- 
sonably common' in their mock galaxy catalogues derived 
from TV-body simulations. Here we take a more cautious ap- 
proach and admit the possibility of a systematic difference 
between the NGP and SGP, possibly a photometric offset. 
Indeed, we have found an offset of 0.056 mag between the 
MGC and the 2dFGRS-IC NGP. Here and in the previous 
section we have included this offset in the conversion from 
bj to Bmgc- Applying this offset only to the NGP, but not 
the SGP, makes the difference worse. 

Since we require homogeneity for the LSS correction 
to work we choose not to use any SGP data in Nic- On the 
other hand, the NGP seems reasonably homogeneous. In the 
top panel of Fig. 1141 we plot as solid dots the counts from 
21 independent, MGC-shaped subregions of the NGP. We 
can see that the MGC region is typical and from Fig. I13l we 
have already seen that the variance among these subregions 
agrees well with that expected. Hence we will use the NGP 
counts as Nic- 

What are the error properties of the MGC-2dFGRS 
counts? Defining f- mc = Nugc/Nic-ugc we have 



2 



2 



Nh + fin 



2 



+ 2pfi nc NlCCTf inc <7N ll 



(11) 



where p is the correlation coefficient between /i nc and Nic- 
A large value of /i nc indicates that many galaxies in the 
MGC region are missed by the 2dFGRS-IC. However, this 
will be true everywhere if the properties of the 2dFGRS-IC 
are homogeneous. Hence Nic will be low if /j nc is large and 
vice versa, i.e. p must be negative and we will err on the side 
of caution by neglecting the last term. The first term can be 
estimated using the jackknife method. For each 'subsample' 
(generated by excluding one field) we calculate N^gc- Since 
./Vic does not vary from subsample to subsample the result- 
ing error estimate is just the first term above, which we plot 
as open diamonds in Fig. 1131 The estimate agrees very well 
with the Poisson errors of Nmgc (or Nuhc)- Since /i, lc ~ 1 
we will use 



2 



Nmgc 
A 



(12) 



where we estimate the second term as well as the rest of the 
covariance matrix using the integral constraint method. The 
result is shown as '+' symbols in Fig. 1131 We find that the 
errors of the MGC-2dFGRS counts are smaller by 12-32 per 
cent than the errors of the counts obtained from the MGC 
alone. 



6.3 Correcting for asteroids 

As a result of their extended morphology aste roids will be 
classified as galaxies in MGC-FAINT. Recently. Ilvezic et al] 
(2001) performed accurate measurements of the shape of 
the number counts of both C- and S-type asteroids for 
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Figure 15. Upper panel: the dots show the number counts of 
all objects in the cosmic ray band delineated by the two solid 
lines in Fig. [TO] The squares show the number counts of objects 
selected by applying the additional cut FWHM < 0.9 X seeing. 
Lower panel: the ratio of the above counts appears to remain 
constant to B MGC = 22.7 mag (dotted line), where real objects 
begin to contaminate the counts from the cosmic ray band. The 
incompleteness of cosmic rays selected by the additional FWHM 
cut is found to be 18.1 per cent and is indicated by the solid 
horizontal line. 

Bmgc ^22.5 mag from SDSS-CD. Here we adopt their 
model, normalized to the number of visually identified as- 
teroids in MGC-BRIGHT, and subtract the model asteroid 
counts from the galaxy counts (cf. Fig. 1161 . 

6.4 Correcting for cosmic ray incompleteness 

In Section 14.5.31 we noted that the identification of cosmic 
rays is incomplete for B M gc — 22.7 mag. Here we will esti- 
mate the size of the incompleteness and correct the galaxy 
number counts accordingly. 

In the upper panel of Fig. 1151 we plot as solid dots the 
number counts of all objects in the cosmic ray band delin- 
eated by the solid lines in Fig. 1101 For Bjj GC < 22.7 mag 
these represent the complete cosmic ray counts but be- 
yond 22.7 mag they are 'contaminated' by real objects. At 
21 < B' MGC < 22 mag the counts are very steep and would 
exceed the galaxy counts by B M gc ~ 23 mag if they did 
not flatten out. In the range 22 < B' MGC < 22.5 mag the 
counts indeed show some flattening. This behaviour makes 
it impossible to reliably extrapolate the cosmic ray counts 
to fainter magnitudes (as we do for the stars and asteroids) . 

We plot as solid squares the number counts of those 
objects that lie in the cosmic ray band and have FWHM 
< 0.9 x seeing, which we use as cosmic ray identification 
criteria at B M gc > 22.7 mag. By taking the ratio of the 
two sets of counts in the lower panel of Fig. I15l we find that 
for Bmgc < 22.7 mag 81.9 per cent of all objects within 
the cosmic ray band also have FWHM < 0.9 x seeing and 



Figure 16. MGC number counts as indicated. The open circles 
refer to the galaxy counts corrected for LSS (£?MGC < 19 ma g)i 
asteroids (Bmgc > 20 mag) and cosmic ray incompleteness 
(-B M GC > 22.7 mag). The open triangles show the cosmic ray 
counts corrected for incom pleteness at -B( Jf -.^ > 22.7 mag. The 
solid line shows the adopted llvezic et alj 1200 ll) model for the as- 
teroid counts, which we extrapolate beyond the faint limit probed 
by these authors (dashed line). 

that this ratio remains constant to within the errors. We 
now assume that this ratio remains constant at -B M gc — 
22.7 mag and correct the cosmic ray number counts in the 
range 22.5 < Bmgc < 24 mag by applying a factor of 1.22. 
The difference between the corrected and uncorrected counts 
is then subtracted from the galaxy counts. 

6.5 Number counts 

We present the MGC number counts in Fig. Ilbl and list the 
galaxy counts in Table [5] Note that the correction for LSS 
(which affects the counts only in the range 16 < Bmgc < 
19 mag) is relatively large at -Bmgc = 17.25 mag (17 per 
cent) but < 4 per cent in all other bins (cf. middle panel of 
Fig. 1131 . At < 1.1 and ~ 1.4 per cent the asteroid and cosmic 
ray incompleteness corrections (which affect the counts only 
at Bmgc > 20 and B' MGC > 22.7 mag, respectively) are 
smaller still. 

Fig. 1171 compares the MGC galaxy number counts to 
those of previous surveys. In the range 16 < Bmgc < 24 mag 
our counts lie among the montage of previous publica- 
tions and provide a fully consistent, uniform, well-selected 
and complete sample spanning eight magnitudes. The data 
thus represent a significant connection between the local 
photographic surveys (B < 20 mag) and the deep pencil 
beam CCD surveys (B > 22 mag). At bright magnitudes our 
counts are consistently higher than the original APM counts 
( Mad dox et al.lll990bl) and suggest that the steep rise of the 
APM counts is an artefact. We note that the 2dFGRS-IC 
is a substantially revised version of the APM catalogue and 
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Table 2. MGC galaxy number counts. 



BmGC 
(mag) 


Af a 


Area 
(deg 2 ) 


MGC 


N 


yycor cL 
MGC 




16.25 


78 


30.84 


2.53 


0.34 


2.44 


0.30 


16.75 


143 


30.84 


4.64 


0.49 


4.90 


0.41 


17.25 


242 


30.84 


7.85 


0.68 


8.67 


0.56 


17.75 


523 


30.84 


17.0 


1.1 


16.77 


0.84 


18.25 


925 


30.84 


30.0 


1.7 


28.7 


1.2 


18.75 


1497 


30.84 


48.5 


2.3 


49.2 


1.6 


19.25 


2433 


30.84 


78.9 


3.2 


78.9 


3.2 


19.75 


4016 


30.84 


130.2 


4.6 


130.2 


1.6 


20.25 


6495 


30.25 


214.7 


6 5 


212.4 


6.5 


20.75 


10568 


30.25 


349.4 


9.2 


346.1 


9.2 


21.25 


17727 


30.25 


586 


13 


581 


13 


21.75 


30074 


30.25 


994 


20 


988 


20 


22.25 


51663 


30.25 


1708 


30 


1698 


30 


22.75 


88353 


30.25 


2921 


45 


2869 


45 


23.25 


129936 


26.71 


4865 


68 


4782 


68 


23.75 


22068 


2.91 


7581 


219 


7451 


219 



"Number of galaxies in the MGC in units of 
^Number counts in units of (0.5 mag) -1 deg 
c Integral constraint estimate of Section l6.ll 
^Number counts corrected for LSS (16 < 
asteroids (-Bmgc ^ 20 mag) and cosmic 
( B MGC ^ 22 - 7 ma §)- 



(0.5 mag) 
-2 



-Bmgc < 19 mag), 
ray incompleteness 



shows a less pronounced rise. Although some local effect is 
still evident (cf. Fig. I18H there is no more need for strong lo- 
cal evolution of the luminous galaxy p opulation as originally 
put forward bv lMaddox et all l)l990br) . 

Fig. 1181 shows a close-up comparison of the MGC- 
B RIG HT, SDSS-EDR, 2dFGRS-IC NGP and lGardner et all 
( 1996) counts, where we have normalized the counts by an 
arbitrary linear model in order to make any differences more 
easily discernible. 

We have used the SDSS-EDR g* counts of the North- 
ern stripe of lYasuda et all l)200lf) and conv erted to -Bmgc 
using equation , the colour equation of iFukueita et all 
( 1996 ) and the galaxy colours as a function magnitude 
ofJ \as_uda_et all a200lf) . Alternatively, we could have used 
lYasuda^tHdT si? counts which were derived using individual 
galaxy colours rather than a mean colour for a given mag- 
nitude. However, the colou r equation used bvjYasuda et all 
is inconsistent with that of IFukueita et alJ il996ft . We pre- 
fer the latter over the former because it was used in a 
direct comparison of SDSS-EDR and MGC magnitudes, 
which showed essentia lly no overall zer o-point offset be- 
tw een the two su rveys JCross et all 120031) . It was also used 
in iNorberg et all s comparison of the 2dFGRS and SDSS- 
EDR photometry, which also showed only a small difference 
of Am(2dFGRS— SDSS-EDR) = 0.058 mag. 

Fig. EH1 shows that the MGC and SDSS-EDR counts 
are in reasonably good agreement. In particular, both show 
a characteristic change of slope near -Bmgc = 18 mag. How- 
ever, the SDSS-EDR counts are higher than the MGC counts 
in t he range 18 < -Bmg c < 21 mag by 4-7 per cent. Given 
that ICross et alJ ~ j2003T) found no offset between the MGC 
and SDSS-EDR photometries this difference cannot be ex- 
plained photometrically. However, the error bars in Fig. 1181 
which include an LSS component, indicate that the differ- 
ence is not significant and may well be caused by LSS. 

Nevertheless, we point out that the definition of a 
galaxy varies among different SDSS publications. The cen- 
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Figure 17. The galaxy number counts derived from the MGC as 
compared with the number counts of various other authors. The 
various counts were converted onto t he -Bmgc system assuming 
(B - V) = 0.94 iNorberg et al . 2002). The 2dFGRS no-evolution 
model counts are calculated from the luminosity function param- 
eters given in Table El (cf. Section ITU. Also shown is the 0.6 
'Euclidean' slope. 
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Figure 18. A close-up comparison over the important normal- 
ization range 16 < -Bmgc < 20 mag of the MGC number counts 
with three other datasets as indicated. The MGC error bars are 
the integral constraint estimates (cf. Sections 16.11 16.21 and Ta- 
ble 0. The SDSS-EDR and 2dF GRS-IC errors are al so integral 
constraint estimates, whereas the lGardner et all Jl996f) errors are 
those published. 
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tral tool of the SDSS-EDR star-galaxy separation is the dif- 
ference between the point spread function (PSF) and model 
magnitudes, c, of an object. The standard SDSS-EDR pro- 
cedure is to classify an object as a galaxy if i t satisfies the 
condition c > 0.145 iStouehton et aljliool ). lYasuda et alJ 
use a slight variant in that they requi re c > 0.145 in two o f 
the three bands g* ,r* and i* . However. lBlantonej i _alj j|200l[) 
use a cut-off value of 0.242 in the r* band~ lstoau^set^rt] 
( 2002) even use a value of 0.3 in the SDSS spectroscopic 
target selection in addition to a surface brightness cut and 
further selection based on photometric flags. They find 
that at the bright magnitudes of the spectroscopic sam- 
ple (r* < 17.8 mag) only 10 per cent of all objects with 
0.15 < c < 0.3 are actually galaxies (see also their fig. 7) 
and they conclude that their sample completeness is > 99 
per cent. 

Hence contamination by misclassified objects may in 
part be responsible for the observed difference between the 
MGC and SD SS-EDR counts. In deed, the object-by-object 
comparison of lCross et alJ J2003L in the range 16 < Bmgc < 
20 mag) found the SDSS-EDR galaxy sample contaminated 
by stars and artefacts at the ~ 1 per cent level. 

Comparing the uncorrected MGC and 2dFGRS-IC 
counts we note that the latter are generally lower by 5- 
10 per cent at Bmgc < 18.5 mag, although the difference 
is again only marginally significant. We note that in the 
tra nsformation from bj to - Bmgc magnitudes we have used 
the iBlair fc Gilmord (Il982l) colour term of 0.28 (cf . equa- 
tion !13> and a global mean galaxy colour of (B — V) = 0.94 
iNorberg et al.ll2002l). Using ins t ead a colour term of 0.35, 
as favoured bv iMetcalfe et all jl995l) . or a larger mean 
galaxy colour, which may be reasonable for the fainter 
galaxies, both exacerbate the difference. In any case, the 
small photometric offset between the MGC and 2dFGRS- 
IC photometries has already been taken into account and 
hence the difference cannot be explained photometrically. 
However, a difference of this magnitude is easily explained 
by the incompleteness and stellar contamination in the 
2dFGRS-IC, whi ch were found to be ~ 9, and ~ 6 per 
cent, respectively JPimbblet et alfeOOltlNorberg et all200l 
ICross et al.ll2003l) . 

As far as we are aware, apart from the SDSS-EDR 
and the MGC the largest CCD survey from which num- 
b er count data have be en published is the 8.5 deg 2 survey 
of lGardner et all dl996T) . Although generally lower than the 
MGC counts by ~ 10 per cent the error bars indicate that 
the difference may well be due to LSS. 



7 REFINING THE FIELD GALAXY 
LUMINOSITY FUNCTION 

The field galaxy luminosity function (LF) ha s been mea- 
sured by many surveys but, as pointed out bv lCross et alJ 
( 2001), their results - and in particular their derived normal- 
izations - are inconsistent with each other (cf. also Fig. 119b ,'). 
Depending on which of the many LFs one selects, one will 
predict markedly different galaxy counts, even locally (see 
Fig. 119b ) . In the following section we use the local galaxy 
counts derived in Section HJ to assess which LFs predict 
counts consistent with the data and to provide stringent 
constraints on the LF normalization. 



7.1 Modeling the number counts 

In Table |3] we list the LF paramete rs from the 
2dFG RS jNorberg et alJ I2002D SPSS -CD (iBlanton et all 
l200ll) . ESP JZucca et al.1 Il997l) CS ferown et al.l l200lF 7 
Purham/UKST jRatcliffe et alJ ll99Sl). Mt St r omlo/ APM 
lLovedav etafl I1992T) SSR S2 llMarzke et all 119981) and 
NOG jMarinoni et al]ll999l) surveys. All M* values have 
been converted to the Bmgc system using equations 1131 
below. These wer e derived from e quati on l |5^ and the 
colour equations ofl Blair fc GilmoTel lll982h .lFuku^h^a_e_t_aT' 



(1996) (SDSS-CD), Bro wn et al.l ll200ll) (TJS). lAlonso et al 



( 1994 1) (SSRS2 ) iKirshner. Oemler fc Schechterl (11978!) and 
iPeterson et alJ jl986h (NOG). We also correct for a known 
zero-point offset in the 2dFGRS photometry: 3 

Bmgc = &2ctfgrs - 0.056 + 0.14(5 - V) 

Bmgc = 6j + 0.14(5 - F) 

Bmgc = g* + 0.12 + 0.30(B - V) 

Bmgc = &SSRS2 - 0.02 - 0.145(B - V) (13) 

Bmgc = Vbs - 0.072(1/ - R) + 0.8553(B - V) 

Bmgc = B RC 3 + 0.74 - 0.51(5 - V) 

We assu me (B -V) = 0.94 llNorberg et alJ 120021) and 
(V-R) = 0.53 jBrown et al.ll200ll) . 

Where available we have used the parameters for the 
currently favoured (JJm.Oa) = (0.3,0.7) cosmology (2dF- 
GRS, SDSS-CD and CS). The N OG parameters do no t de- 
pend on the cosmological model (iMarinoni et al]ll998l) . the 
SSRS2 used (0.4,0) and the rest used (1,0). We have at- 
tempted to correct the M* values of these latter surveys to 
a (0.3, 0.7) cosmology using 



M* (0.3, 0.7) = M*(n M ,0 A ) +5 log 



r L (^M, ^a; z) 
r L (0.3,0.7; z) 



(14) 



where tl and z are the luminosity distance and median red- 
shift of the survey, respect ively. Since all of these surveys 
fixed (j>* by using one of the lDavis fc Huchral l)l982f) estima- 
tors of the mean galaxy density we also corrected their <f>* 
values using 



0* (0.3, 0.7) = <j}*(n M ,n A ) 



^j(S7m, ^a; z) 
% (0.3, 0.7; z) 



(15) 



where dV/dz is the comoving volume element. Since the 
2dFGRS, SDSS-CD and CS have all derived their LF pa- 
rameters for several different cosmological models we can 
use these to test this simple correction procedure. In all 
cases we find good agreement (well within the quoted error 
bars) between the transformed and measured values. The 
final corrected Bmgc LFs are shown in Fig. I19f a') . 

3 Although the MGC was used during the 2dFGRS photom- 
etry recalibration procedure (100k public release) to establish 
non-linearity corrections, it was never used in the determina- 
tion of the overall zero-point of the 2dFGRS. Hence a zero- 
point offset between the two surveys is possible. We have 
performed an object-by-object comparison of the photometry 
of 5996 galaxies in common to the two surveys and found 
A m(MGC-2dF G RS) = (-0.056±0.005) mag. Similarly, combin- 
ing Eros^t^ijsj2003) Am(MGC— SDSS-EDR) = -0.002 mag 
with lNorberg et all s Am(SDSS-EDR— 2dFGRS) = -0.058 mag 
gives Am(MGC-2dFGRS) = -0.060 mag. 
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Table 3. Parameters of published luminosity functions (cf. Fig. 1191 and their revised normalizations and luminosity densities. For each 
LF we also list the probability that its number count predictions fit the data before and after the adjustment of <j>* (eight and seven 
degrees of freedom, respectively). 
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! ) 
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10 


~ 2 h 3 Mpc" 








(10 s h 


L 


Mpc- 3 ) 


2dFGRS 


i>2dFGRS + 0-07 


0.1 


-19.59 


± 


0.008" 


-1.21 


± 


0.01 


1.61 


± 


0.06 


0.18 


1.69 


± 


0.04 ± 0.03 


± 


0.03 


0.44 


2.03 


± 


0.05 


SDSS-CD 


9* + 0.40 


0.1 


-19.64 


± 


0.04 


-1.26 


± 


0.05 


2.06 


± 


0.23 





1.56 


± 


0.04 ± 0.11 


± 


0.03 


0.40 


2.06 


± 


0.11 


ESP 


bj + 0.13 


0.1 


-19.59 


± 


0.08 


-1.22 


± 


0.07 


1.65 


± 


0.3 


0.45 


1.68 


± 


0.04 ± 0.20 


± 


0.03 


0.42 


2.04 


± 


0.12 


cs 


Vas + 0.77 


0.064 


-19.53 


± 


0.09 


-1.09 


± 


0.09 


1.87 


± 


0.21 


0.62 


1.91 


± 


0.04 ± 0.19 


± 


0.06 


0.60 


1.96 


± 


0.10 


Dur./UKST 


bj + 0.13 


0.052 


-19.61 


± 


0.10 


-1.04 


± 


0.08 


1.53 


± 


0.3 


9 X 10~ 7 


1.78 


± 


0.04 ± 0.23 


± 


0.06 


0.77 


1.90 


± 


0.09 


Str./APM 


bj + 0.13 


0.051 


-19.43 


± 


0.13 


-0.97 


± 


0.15 


1.26 


± 


0.15 





2.20 


± 


0.05 ± 0.37 


± 


0.07 


0.67 


1.90 


± 


0.13 


SSRS2 


&SSRS2 - 0.16 


0.02 


-19.61 


± 


0.06 


-1.12 


± 


0.05 


1.24 


± 


0.2 





1.73 


± 


0.04 ± 0.14 


± 


0.07 


0.65 


1.95 


± 


0.09 


NOG 


Srcs + 0-26 


0.01 


-19.80 


± 


0.11 


-1.11 


± 


0.07 


1.40 


± 


0.5 


0.87 


1.42 


± 


0.03 ± 0.20 


± 


0.08 


0.82 


1.89 


± 


0.14 



a The three quoted errors are due to: (i) statistics and LSS, (ii) the correlated errors on M* and a and (iii) the uncertainty in k + e- 
corrections. 



6 We have reduced INorberg et all s zero-point uncertainty from 0.04 to 0.005 mag because we have found that the 2dFGRS and MGC 
photometries agree to this level once the zero-point offset has been applied. We have also excluded the error induced by k + e-corrections 
as we treat this error separately and consistently for all surveys in Section l7.3l 



To model the counts we use (Am, ^a) = (0.3, 0.7), a k- 
correction of 2.5z and an e-correction of 2.5 log[(l ± z) -0,75 ]. 
At z < 0. 4 this combinat i on ma tches well the fe+e-correction 
given bv INorb erg et all j2002|) f their fig. 8) which was de- 
rived using iBrtizued^^Cnarion 1119931) models to match the 
colour-redshift trend seen in the 2dFGRS (using SDSS-EDR 
colours). In principle, for each of the surveys in Tabic [jj|wc 
should use the same fc(±e)-correction as was used in the 
derivation of the Schechter parameters of that survey or else 
correct for the use of a different k + e. However, in several 
cases different fc-corrections were used for different galaxy 
types, sometimes interpolating between types, and hence a 
mean correction is not readily available. 

Fig. I19f b) shows how the counts predicted by the vari- 
ous LFs compare with the MGC-BRIGHT counts (corrected 
for LSS). Clearly, the data prefer some models over others. 
We have performed a goodness-of-fit test for the predicted 
counts in the range 16 < Bmgc < 20 mag (using the full 
covariance matrix of the counts) and list the resulting prob- 
abilities in column 7 of Table [3] The model counts based on 
the Durham/UKST, Mt Stromlo/APM and SSRS2 LFs are 
all significantly too low, while the SDSS-CD model counts 
are too high. The other models provide good fits. These 
conclusions do not change if we compare the models to the 
uncorrected counts. 

From Fig. I19f b) it is clear that the main difference be- 
tween the various predicted counts and the data lies in their 
normalization, Nis, and not their shape. As discussed in 
Section it is precisely this uncertainty in Nis that is com- 
monly referred to as the 'normalization problem'. 



7.2 The normalization problem 

For a given set of Schechter parameters what is the uncer- 
tainty in Nig induced by the error of each of the parameters? 
In the general case there is no analytic formula relating iVis 
with the Schechter parameters. Hence the 'Euclidean' case 
must serve as a guideline. Ignoring binning effects we have 
Ai8 oc 4>* L* 3/,2 r(a + |). Using the intersurvey rms of each 
parameter as an indication of its uncertainty we find 



dN 1: 



■A<j>* 



dN 18 
dM* 



AM* 



dN ls 
da 



Aa = 13 : 10 : 1 



(16) 



and similar or worse values for the individual surveys. Hence 
we identify the errors in both (f>* and M* as the main causes 



of the uncertainty in iVi8. What are the sources of these 
uncertainties? 

7.2.1 Photometric error 

In addition to the difficulty of reliably cali brating the 
photo metry of large photographic surveys ICross fc Driverl 
(2002) showed that the derived LF parameters depend 
on the limiting isophote and the photometric method 
(e.g. isophotal, corrected, total, etc.). They concluded that 
when using isophotal magnitudes for a limiting isophote of 
25 mag arcsec -2 one might expect errors in M* of up to 
±0.4 mag, in <j>* of up to ±10 per cent and in a of up to 
±0.01. Hence photometric uncertainties contribute substan- 
tially to the normalization problem. 

7.2.2 Incompleteness 

Although the underestimation of the magnitudes due to the 
limiting isophote will cause some galaxies to fall below the 
limiting magnitude of a survey, it also causes the derived 
volume over which such galaxies can be seen to be under- 
estimated. In addition, Fig. 1121 (together with visibility the- 
ory) shows that at Bmgc ^ 19 mag apparently no significant 
galaxy population below 25 mag arcsec -2 exists. In other 
words, at these surface brightness limits it is typically a case 
of missing li ght from the outer isophotes rather than miss- 
ing galaxies: ICross et al.l i2003l) found only ~ 0.5 per cent of 
galaxies missing from the 2dFGRS due to their low surface 
brightness. However, blended and unresolved objects take 
the total incompleteness of the 2dFGRS to ~ 9 per cent, 
which may be indicative of photographic surveys in general. 

7.2.3 Large-scale structure 

LSS will clearly aff ect all but the very largest surveys. 
INorberg et alJ i2002f) found that the 173 deg 2 overlap region 
between the SDSS-EDR and the 2dFGRS was overdense by 
5 per cent relative to the full 2dFGRS. Similarly they found 
a 7 per cent difference (to fcj = 19.2) between their 740 deg 2 
NGP and 1095 deg 2 SGP regions. Hence it is clear that much 
shallower surveys like the SSRS2 (&n m = 15.5 mag) cannot 
accurately measure <f>* even if they have a very large angu- 
lar size. On the other hand, due to the filamentary structure 
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Figure 19. The left-hand panels show local LFs before (upper) and after (lower) renormalization to the corrected MGC-BRIGHT 
counts. The right-hand panels show the number count predictions of these LFs and the LSS-corrected MGC-BRIGHT counts (i.e. the 
MGC-2dFGRS counts at Bmgc < 19 mag). After the recalculation of the LF normalizations, <j>* , we note that all LFs match the data 
equally well. 



of the Universe even deep surveys are susceptible to LSS if 
they extend only over a small solid angle. 



suits iYasuda et al .l200lHNorberg et al l2002HBlanton et alJ 

2003). 



7.2.4 k(+e) -corrections 

Finally, we remark that each of the surveys listed in Table 
louses different fc-corrections, and only the 2dFGRS use e- 
cor rections. Indeed, evolu tionary effects, which were ignored 
by iBlanton et all (1200 ll) . are the reason why the SDSS- 
CD normalization is so high. When evolution is included 
(or when normalizin g to the S DSS-ED R co unts instead of 
using the method of Davis fc Huchralll982T) the SDSS-CD 
normalization agrees much better with the 2dFGRS re- 



7.3 Determining (f>* 

Although both M* and cj>* appear to contribute equally to 
the uncertainty in Ni$ it is customary (and sensible) to fix 
M* and use Nig to constrain (f>* alone instead of a combi- 
nation like 0*Z/* 3 / 2 . Over the range 16 < Bmgc < 20 mag 
the MGC arguably provides the most reliable number count 
data in existence and hence it provides very reliable con- 
straints on <f>* . We now determine 4>mgc as a function of M * 
and a by fitting the model counts for a given combination of 
M* and a to the corrected MGC-BRIGHT counts over the 
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Figure 20. The contours show the best-fitting value of the LF 
normalization, <Pmgc ^ or gi ven values of M* and a. These were 
derived by x 2 -minimization over the range 16 < Bmgc < 20 mag 
(cf. Fig. 1191 using the full covariance matrix). The values of the 
contours increase from 0.011 h 3 Mpc - 3 on the left in steps of 
0.001 to 0.027 on the right as indicated. The error ellipses of the 
LFs of Tablc[3]in the M*-a plane are also indicated (assuming a 
correlation coefficient of 0.75). The corresponding best-fit </> m gc 
values are listed in Table [3] The light and dark shading in the 
upper right corner marks those regions in the M*-a plane where 
the shape of the predicted counts disagrees with the data at the 
90 and 95 per cent levels, respectively. 



range 16 < Bmgc < 20 mag (i.e. the MGC-2dFGRS counts 
at Bmgc < 19 mag), where </>* is the free parameter. We 
use the full covariance matrix for the fits. Fig. I20l shows the 
result as a contour plot in the M*-a plane. 

In particular, we derive the </>mgc values appropriate 
for the LFs listed in Table [3] The revised values and their 
associated probabilities are listed in columns 7 and 8 of that 
table. Figs. I19f c'l and (d) show the resulting LFs and their 
predicted counts after revision. All surveys now give compa- 
rable and reasonable probabilities. Note, however, that the 
variation of <j>* among the different surveys after renormal- 
ization is no less than before. This is due to the fact that 
the counts constrain the combination 0*L* 3,/2 and that the 
variation in £* 3 / 2 i s comparable to that of 4>* before renor- 
malization (as seen from equation 1161 . 

In Table |H] we also list the errors on <^> M gc induced by 
various sources. First, we list the statistical error derived 
from the x 2 fit. This error includes remaining uncertainties 
due to LSS because we have used the full covariance matrix 
of the counts in the fit. Secondly, we list the errors due to 
the uncertainties in M* and a, which can be essentially read 
off from Fig. 1201 Here we take into account that the errors 
on M* and a are generally hi ghly correlated by as suming a 
correlation coefficient of 0.75 (Bl anton et al-llioOlT) . Finally, 
we have estimated the errors resulting from uncertainties 
in the k + e-correction. This error is slightly more subtle 
than the previous ones. For example, increasing k + e has 
the effect of moving the contours in Fig. 1201 to the left and 
up but it also moves the LFs in a similar direction, thus 
again decreasing the origin al change in 0mgc f° r a given 
survey. INorberg et al] i2002T) estimated errors on their k + e- 
correction by demanding statistical consistency between the 



LFs derived from their high and low-z samples. Since our 
k+e is essentially identical to theirs we also adopt their error 
of ±18 per cent. We estimate the resulting <^> M gc error for a 
given survey by first adding ±Q.18(fc + e)(z) to its M* value 
and then re-deriving (/>mgc while using the increased (or 
decreased) k + e. As in Section lV.ll we thus approximate the 
effect on M* by evaluating A( k + e) at the me dian redshift of 
the survey. Comparison with INorberg et alt s k + e induced 
error on M* indicates that this procedure underestimates 
the effect on M* , and hence we overestimate the error on 
^mgc ■ We also conservatively ignore any effect on a, which 
is difficult to gauge. 

7.4 The local luminosity density 

Taking M@ b = +5.3 mag and using the revised LF nor- 
malizations we calculate the luminosity density, ji 3 = 
4>*Ll 3 T{a + 2). Note that we convert MJ MGC to M 6 * r The 
calculated values of jbj are listed in the last column of Ta- 
ble []] where the errors include all the M* , a and </>mgc 
uncertainties while taking into account the correlation be- 
tween the errors of M* and a (with a correlation coefficient 
0.75 as above) as well as the virtually perfect correlations 
between the second error component of 0mgc an< l l ne er- 
rors of M* and a and between the k + e errors on M* 
and 0mgc- We find a weighted mean 6j luminosity den- 
sity of = (1.986 ± 0.031) x 10 s h L Mpc -3 . Note 
that the 2dFGRS and SDSS-CD, which must be consid- 
ered the most reliable LFs, give a slightly higher value of 
j b] = (2.035 ± 0.046) x 10 8 h L e Mpc -3 . 

7.5 Constraints on M* and a 

In the previous section we minimized x 2 as a function of <jf 
only while treating M* and a as fixed parameters. (f>* only 
affects the normalization of the predicted counts but not 
their shape which, in principle, also contains information. 
However, for the LFs listed in Table |3 the fits are already 
statistically acceptable (cf. column 9 of Table |3J and hence 
we cannot expect to derive useful limits on any other pa- 
rameters. 

Nevertheless, in Fig. 1201 we mark as grey-shaded regions 
those combinations of M* and a where the adjustment of 
(f>* does not result in an acceptable fit (at the 90 and 95 
per cent confidence levels), i.e. where the shape of the pre- 
dicted counts disagrees with the data. As expected these 
limits cannot exclude any of the measured values. 



8 CONCLUSIONS 

Here we have presented a detailed description of the 
Millennium Galaxy Catalogue (MGC), a deep (nu m — 
26 mag arcsec -2 ) wide (37.5 deg 2 ) survey along the equa- 
torial strip from 9 h 58 m to 14 h 47 m . We have demonstrated 
that the internal photometric accuracy of the MGC is 
±0.023 mag and that the astrometric accuracy is ±0.08 arc- 
sec in both RA and Dec. Using SExtractor we have de- 
rived a source catalogue containing over 1 million objects 
spanning the range 16 < -Bmgc < 24 mag. All non-stellar 
detections brighter than Bmgc = 20 mag have been visu- 
ally inspected and the objects repaired where necessary. We 
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have taken care to exclude objects from regions where the 
photometry is likely to be erroneous, resulting in a robust 
and clean estimation of the galaxy number counts over the 
range 16 < -Bmgc < 24 mag. These data finally connect the 
faint pencil beam CCD surveys of the past decade to the 
local Universe. The selection boundaries of the MGC are 
well defined and to Bmgc = 20 mag we are demonstrably 
robust to star-galaxy separation and low- and high-surface 
brightness concerns. We contest that the MGC galaxy num- 
ber counts in this range are the state of the art, superseding 
all previous intermediate number count data. 

We use the counts to test various estimates of the galaxy 
luminosity function and find that many of them predict 
counts where the normalizations are inconsistent with our 
observations. In Fig. !2Ul we present the best-fitting value of 
</>* as a function of M* and a. In Table0we list the appro- 
priate cj>* values for a number of popular B-band luminosity 
function estimates. 

From these revised values we constrain the &j-band lu- 
minosity density of the local Universe for each of these 
luminosity functions. We find jtj = (1.986 ±0.031) x 
10 8 h L Q Mpc -3 . The 2dFGRS and SDSS-EDR consis- 
tently give a slightly higher value of j^j = (2.035 ± 0.046) x 
10 s hL @ Mpc" 3 . 
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